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ABSTRACT 

Zeeman observations of molecular clouds yield the line-of-sight component 
^ Blos of the magnetic vector B, which makes it possible to test the two major 

>- ■ extreme-case theories of what drives star formation - ambipolar diffusion or tur- 

^ ■ bulence. However, only one of the three components of B is measurable, so tests 

OO ' have been statistical rather than direct, and they have not been definitive. We 

report here observations of the Zeeman effect in the 18-cm lines of OH in the 
(^ ■ envelope regions surrounding four molecular cloud cores toward which detections 

52 ' of Blos have been achieved in the same lines, and evaluate the ratio of mass to 

magnetic flux, M/$, between the cloud core and envelope. This relative M/$ 
1^ ■ measurement reduces uncertainties in previous studies, such as the angle between 

j_| ' B and the line of sight and the value of [OH/H]. Our result is that for all four 

clouds, the ratios 7?. of the core to the envelope values of M/<l> are less than 1. 
Stated another way, the ratios TZ' of the core to the total cloud M/$ are less than 
1. The extreme case or idealized (no turbulence) ambipolar diffusion theory of 
core formation requires the ratio of the central to total M/$ to be approximately 
equal to the inverse of the original subcritical M/$, ot TZ' > 1. The probability 
that all four of our clouds have 7?.' > 1 is 3 x 10~^; our results are therefore sig- 
nificantly in contradiction with the hypothesis that these four cores were formed 
by ambipolar diffuson. Highly super- Alfvenic turbulent simulations yield a wide 
range of relative M/$, but favor a ratio 7?. < 1, as we observe. Our experi- 
ment is limited to four clouds, and we can only directly test the predictions of 
the extreme-case "idealized" models of ambipolar-diffusion driven star formation 
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that have a regular magnetic field morphology. Nonetheless, our experimental 
results are not consistent with the "idealized" strong field, ambipolar diffusion 
theory of star formation. Comparisons of our results with more realistic models 
and simulations that include both ambipolar diffusion and turbulence may help 
to refine our understanding of the relative importance of magnetic fields and 
turbulence in the star formation process. 

Subject headings: ISM: clouds, evolution, magnetic fields — polarization — stars: 
formation 



1. Introduction 



Understanding star formation is a fundamental astrophysical problem. iMcKee fc Ostriker 



( 120071 ) have comprehensively reviewed the field. For thirty years what has sometimes been 
called the "standard model" has been that magnetic fields control the formation and evolu- 
tion of the molecular clouds from which stars form, including the formation of cores and their 
gravitational collapse to form protostars. However, in recent years doubts about the validity 
of this model have been raised by those who argue that turbulence controls the formation of 
clouds and cores, with cores either dissipating back into the general interstellar medium or 
collapsing and forming stars if they are self-gravitating when formed. In spite of decades of 
intense research, there is still not consensus on the role that magnetic fields play in the star 
formation process. 

Detailed theoretical w ork on the strong mag netic field "standard model" h as been carried 



out by a number of groups: IShu et al.l (119991 ) and lMouschovias fc Ciolekl (119991 ) have reviewed 
and summarized the state of this theory. The fundamental principle is that clouds are 
formed with subcritical masses (M < M$ = $/27rvG). Here $ is t he magnetic fiux and 



G is t he gravitational constant, and the expression for M$ is from iNakano fc Nakamura 



(|l978f): other expressions for M$ differ slightly from this depending on cloud structure (e.g.. 



Mouschovias fc Spitzerlll976l ). The magnetic field is frozen only into the ionized gas and dust; 
neutral gas and dust contract gravitationally through the field and the ions, increasing mass 
(but not to first order fiux) in the cloud cores. This process is known as ambipolar diffusion. 
When the core mass reaches and exceeds M$, the core becomes supercritical (M > M$), 
collapses, and forms stars. The magnetic flux remains behind in the envelope. Because the 
ambipolar diffusion time scale for the formation of supercritical cores is fairly long (> 10^ 
yr), molecular clouds would have long lifetimes. The star formation efficiency (the ratio of 
mass in stars in a molecular cloud complex to the interstellar mass) is low (as observed) due 
to the slow rate of star formation and to the fact that much of the mass of the molecular 
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cloud is left behind in the subcritical envelope as the core collapses. 

The idea that star formation is primarily regulated by ambipolar diffusion was standard 
for many years. However, doubts about the validity of this ass umption were raised by th e 
development of a weak- field, super- Alfvenic model of dark clouds (jPadoan fc Nordlundlll999l ) . 
The new weak-field theory had molecular clouds being intermittent phenomena, with short 
(~ 10^ yr) lifetimes. In this theory clouds form at the intersection of turbulent supersonic 
flows in the interstellar medium. Generally, clouds do not become gravitationally bound and 



dissip ate; those that are self-gravitating form stars in essentially a free-fall time (JElmegreen 



2000). Star formation occurs only in the small fraction of the molecular gas that is sufficiently 
dense to be self-gravitating. The star formation efficiency is low due to a small fraction of the 
mass of clouds becoming gravitationally bound. Magnetic fields are present in this theory, 
but they are too weak to be energetically dominant. The rol e of turbulenc e in the energetics 
of the interstellar medium has been a very active area. Padoan et al.l (1200411 pres ented 
evidence in favor of the weak field, super Alfvenic model. lElmegreen fc Scale (2004) have 



written an excellent review of interstellar turbulence, and iMacLow fc KlessenI (120041 ) have 
extensively reviewed arguments that supersonic turbulence controls star formation. 

Early work focused on extreme-case models, either strong magnetic field models that 
did not (at least directly) include the effects of turbulence, or strong turbulence models 
that neglected the effects of ambipolar diffusion. More recent theoretical work has intro- 
duced supersonic turbulence into numerical sta r formation models (e.g. 



Nakamura &: Li 



20051 : iKudoh &: Basull2008l : iNakamura fc Lill2008l ). These more recent models apply to mag- 
netically subcritical regions, and they include the effects of ambipolar diffusion. The role 
of turbulence proves to be important in these models, significantly shortening, for example, 
the ambipolar diffusion timescale. In this paper we focus on a direct comparison of obser- 
vational results with the strong field, ambipolar diffusion model without turbulence. Such 
models have a smooth morphology of the magnetic field which deforms into an hourglass 
shape with axis along the mean field direction, and make an analytic prediction that can be 
tested directly. The more complicated models that come from simulations that include both 
magnetic fields and turbulence can also be tested against the results of this paper, and we 
discuss how such tests might be carried out. 

Our approach to testing star formation theory has been to measure magnetic field 
strengths in molecular clouds in order to see whether they are weak or strong. The cru- 
cial parameter is the ratio of the mass to the magnetic fiux, M/$, which is of course closely 
related to M$. If M/$ is observed to be significantly supercritical, particularly at lower 
densities, the magnetic support theory is not viable. If it is observed to be subcritical at 
lower densities, magnetic fields would be too strong for the intermittent, turbulent theory to 
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hold. The M/$ parameter provides in principle a straightforward, direct test to discriminate 
between the two extreme theories of star formation. 

The Zeeman effect provides the only known method to directly measure magnetic field 
strengths in dense gas. In this paper we first briefly review Zeeman results and argue that 
they have not yet provided an unambiguous discrimination between the two star formation 
theories. We then describe a new experiment to attempt to overcome the limitations of 
existing Zeeman results and to test directly a prediction of the "standard model" . Finally, 
we describe the observations made for this experiment and the results, and discuss the 
implications of those results. 



2. The Zeeman Effect 

2.1. The Technique 

If a spectral line forming region is permeated by a magnetic field B, the line is split by 
the normal Zeeman effect into three separate frequencies, z/q — z/^, uq, and z/q + z/^, where 
Uz = B X Z , B is the magnitude of B, and Z is the Zeeman coefficient in Hz//iG (note that Z 
is often defined as two times our Z). Unfortunately, most molecules do not have large values 
of Z. In general, only those molecules with an unpaired outer electron will have a Z of order 
the Bohr magneton, Mb = eh/Anmc = 1.40 Hz//iG. For the 1420 MHz H I line, Z = lA 
Hz//xG, so if i? = 1 /iC, the total splitting of the two circularly polarized components would 
be 2.8 Hz, or 6 x 10~^ km s~^. This very small magnitude of the Zeeman splitting makes 
Zeeman observations difficult and consuming of large amounts of telescope time. For most 
molecules, Z will be of order the nuclear magneton, which is 1840 times smaller than Mb- 
Except for the very strong H2O masers, all Zeeman work has involved the few species with 
large Zs: H I, OH, CN, CH, CCS, and SO, with Zeeman detections to date in only the first 
three species. 

The observed data are Stokes parameter I and V spectra for each position. From these 
data we infer the column density N and the line-of-sight component B^os of the magnetic 
vector B. The I and V spectra have channel-to-channel noise that has a Gaussian probability 
density function (PDF) - noise dominated by the stochastic receiver noise. Hence, the 
random error in N is Gaussian, or very nearly so. For the usual case of Zeema n splitting much 



small er than the line width (z/^ << (5z/), only Blos can be determined (e.g., ICrutcher et al. 



19931 ). The Stokes V spectrum actually consists of a possible scaled-down Stokes I signal 
(due to a possible gain difference between the two receivers) and a Zeeman signal that is 
proportional to the first derivative of the Stokes I spectrum. To infer Blos we do a hnear 
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least-squares fit of a x / + 6 x dl/dv to tlie observed Stokes V spectrum, wliere tlie Stokes 
I and V spectra are functions of radial velocity v, and dl/dv is obtained by numerically 
differentiating the Stokes I spectrum. Both of the parameters a and b may be positive or 
negative. Hence, the position and shape of the Zeeman signal are set entirely by the Stokes I 
spectrum; only the amplitude (positive or negative) of the Zeeman signal is a free parameter. 
The sign of b indicates whether the line-of-sight component of the magnetic vector B points 
toward or away from the observer, and the magnitude of b is proportional to the magnitude of 
Blos- The parameters a and b and the mean errors in each parameter come in the standard 
way from this straightforward linear least-squares fitting procedure. Hence, the error in 
Blos that comes from the linear least-squares fit is expected to be Gaussian normal, or very 
nearly so. During the beginning years of Crutcher's observations of the Zeeman effect in OH 
lines, he confirmed by Monte Carlo tests that the linear least-squares fitting procedure gives 
correct results with very nearly Gaussian PDFs. 

Since the magnitude of only one of the three components of the vector B can be mea- 
sured, in general it is necessary to apply statistical techniques in order to determine astro- 
physically meaningful parameters. If one observes a large sample of clouds distributed over 
the sky, the direction of the magnetic fields should be random with respect to the lines of 
sight from the observer. The usual approach is to assume that the total field strength B is 
the same in an observed sample of clouds in order to infer B from measurements of Blos- 



As discussed by iHeiles &: Crutcherl ( 120051 ). for this case both the median and mean values 



of Blos = -B/2, so the mean and rn edian of B may be obtained for the observed sample of 



clouds. As iHeiles fc Crutcherl (120051 ) noted, for PDFs for B other than a delta function, the 
mean and median of Blos ~ -B/2, so useful information about magnetic field strengths from 
Zeeman observations can be obtained even if the PDF of B is unknown ("as is of course the 



case 



Since the magnetic flux is just the magnetic field strength times the spatial area over 
which it is measured, <l> = i? x Area. The mass within this area can be inferred from the 
hydrogen column density, Nh = N {H I) + 2N {H2); allowing for 10% He, Mobs = lAmuNH x 
Area. Thus, M/^ (xNh/B. 



2.2. Previous Zeeman Work 



Most of the earlier Zeeman detections in molecular clouds (e.g., ICrutcherl Il999l ) have 
been toward clouds associated with H II regions. Dark clouds offer the possibility of mea- 
suring the role of magnetic fields at an earlier stage of the star formation process, especially 
for the low-mass star formation case where the "standard model" may best apply. However, 
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until recently there were on ly two dark cloud molecular Zeeman detections (JGoodman et al. 
19891 : ICrutcher et aljEoOsh . 



In order to improve our knowledge of magnetic field strengths in dark cloud cores, 
Troland fc Crutchen (120081 ) used the Arecibo telescope to carry out an extensive program to 
observe the Zeeman effect in the 1665 and 1667 MHz lines of OH. Thirty-three dark cloud 
core positions were observed. They achieved nine detections of Blqs and sensitive upper 
limits for the other positions. Also needed to compute M/$ is an estimate for the column 
density of H2. They obtained this estimate from the OH lines t hemselves. The Arecibo OH 
spectra yield NfOHj . With OH/H = 4 x 10"^ JCrutcheiiEgTgJ ). one can infer N(H2). The 
Troland fc Crutcherl (120081 ) Arecibo maps of OH emission around the core positions showed 
that OH does peak up on the CO and/or NH3 cores, and that OH samples densities up to 
around n(H2) ~ 2 x 10^ cm"^. 



Although the iTroland fc Crutcheii (120081 ) analysis of their data gave the mean value of 
B for their dark-cloud sample, it provided no information about the possible variation of B 
from core to core. They inferred A, the ratio of the observed mean M/<l> to the critical value. 
The data were consistent with the prediction of the strong magnetic field theory - a slightly 
supercritical M/$ if the core morphology is that of a disk, as predicted by the strong-field 
theory. 

The Arecibo Zeeman observations had the potential to eliminate one of the two extreme- 
case theories for the star formation process. If M/$ had been found to be unambiguously 
highly supercritical, the ambipolar diffusion driven theory would have been eliminated. If 
M/$ had been found to be unambiguously subcritical, the turbulence driven theory would 
have been eliminated. The statistical result was that the mean M/$ in molecular clouds was 
observed to be approximately critical, i.e., A ~ 2. Because of the many non-detections, and 
since the OH/H2 ratio is uncertain by a factor ~ 2, a mean M/$ ranging from critical to 
supercritical by ~ 4 was consistent with the data. Of course, there could be a real variation 
with an unknown range in M/$ from cloud to cloud. Hence, although the observations 
have shown that magnetic fields are sufficiently strong that they cannot be ignored, the very 
hard-won observational results cannot rule out either extreme-case theory of star formation. 



3. This Experiment 

The goal of the experiment described in this paper is to perform a definitive test of the 
extreme-case (no turbulence) ambipolar diffusion theory of star formation that circumvents 
the ambiguities inherent in previous Zeeman observational tests. In addition, the experiment 
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will supply valuable observational data against which more complicated models (involving 
both ambipolar diffusion and turbulence) may be tested. This definitive test is measurement 
of the change in M/$ between the envelope and the core of a cloud. The magnetic sup- 
port/ambipolar diffusion theory makes a specific prediction that can be tested. It requires 
that M/$ of the original cloud be subcritical, with ambipolar diffusion accumulating mass 
but not flux in the cloud center, building up a higher density core. Eventually, the core 
becomes supercritical and starts a collapse which is slower than free fall due to the magnetic 
pressure. During this supercritical collapse phase, magnetic flux is dragged inward, so the 
rate of increase of M/$ is slowed. The result is a prediction that the ratio of the core to the 
original cloud M/$ > 1, by approximately the inverse of the amount by which the original 
cloud was subcritical. 

Is the prediction of the ambipolar diffusion theory sufficiently different from the expec- 
tations of the turbulence-driven theory? The answer is clearly yes, for the turbulence theory 
generally predicts the opposite behavior of the mass-to -flux ratio between envelope an d core. 



Simulations of the formation of cores by turbulence (j Vazquez- Semadeni et al.l l2005l ) made 
the usual assumptions for this theory of initially uniform density and magnetic field, driven 
turbulence, and ideal MHD (strict flux freezing, i.e., no ambipolar diffusion). They found 
that M/$ usually decreased with increasing density, the opposite of the ambipolar diffusion 
result. A change in M/$ with ideal MHD seems impossible, but it occurs due to the way 
M/$ is measured. No change in M/$ occurs when an entire flux tube is considered. But 
converging flows perpendicular to B increase B and M within a fixed volume, while flows 
parallel to B can increase or decrease M within that fixed volume. Hence, when measuring 
M/$ within the volume chosen to define a cloud o r core, M/$ can be either larger or smaller 



than the original value. IVazquez-Semadeni et al.l ( 20051 ) noted that for mass and magnetic 



flux conservation, a clump or core within a cloud has a larger density than the mean, but 
a mass smaller than its parent cloud, so a core would generally have a smaller M/$ than 
the parent cloud (since the magnetic flux would be unchanged by core formation but the 
core mass would be only a fraction of the mass within the flux tube). Hence, this turbulent 
simulation result may be general, although of course it would require more studies to verify 
that conclusion. 

Strong magnetic field cases can obtain a similar result. If the mass within a flux tube 
fragments into several cores on a shorter time scale than the ambipolar diffusion one, each 
core will have only a fraction of the mass in the flux tube, an d the M/$ of these cores would 



have decreased with respect to that of the original material (JMouschoviasI Il99ll ) . Although 
the physics driving this fragmentation is different from that considered by the turbulence 
theory, the fundamental reason for the decrease in M/$ is the same, and it is not related to 
ambipolar diffusion. Here we address only testing the ambipolar diffusion driven formation 



of dense cores, for which M/$ must increase. 

A quantitative test of a theoretical prediction requires the design of an experiment to 
test specific numerical models. We designed our experiment to test the "idealized" models 
that have been published by the ambipolar diffusion theorists. We call these models ideal- 
ized because they do not follow the evolution in a consistent way through the formation of 
clouds from a diffuse interstellar medium to the formation of cores in those clouds driven by 
ambipolar diffusion. Instead, the models generally start with an isolated, uniform, spherical 
cloud that is allowed to relax to an equilibrium state, which then evolves due to ambipo- 
lar diffusion. The non-thermal motions and irregular structures of observed clouds, which 
themselves are embedded in a complex interstellar medium, are not directly considered in 
these idealized models. So the models are two-dimensional (with azimuthal symmetry about 
the magnetic field direction) with regular (non-twisted) magnetic fields. To design the ex- 
periment, we looke d at three such idealized amb ipolar diffusion models: (1) a dimensionless 



parameter model (jCiolek &: Mouschoviad Il994j ) for which they listed specific physical pa- 



rameters for comparison with actual molecular clouds; this model with their listed physical 
parameters had an un evolved cloud radius of 4.3 pc; (2) one specifically computed for L1544 
(ICiolek fc Basul 120001) with an une volved cloud radius of about 2.5 pc; and (3) a specific 



model for Bl (ICrutcher et al.lll994l ) with an unevolved cloud radius of 2.9 pc. These radii 
would become the radii of the "envelopes" surrounding the cores formed by ambipolar diffu- 
sion, since the region of the cloud outside the cores would be "held in place" (an ambipolar 
diffusion theorist phrase) by their subcritical magnetic fields. Moreover, the radii of the cores 
in these models were all ~ 0.1 pc. So we need to sample the core with a filled beam of radius 
~ 0.1 pc and the envelope with a beam radius less than ~ 2 pc that excluded the core. A 
second aspect of the idealized ambipolar diffusion models was that the magnetic fields were 
smooth and regular within their unevolved cloud radii, although with an "hourglass" mor- 
phology strongest in the core. Although this regularity of the field was by construction in 
the models, the requirement that M/$ in the clouds be subcritical meant that the magnetic 
fields must be strong, so that even if turbulence had been included in the models, magnetic 
energy would likely dominate turbulent energy and field lines would in fact be quite regular. 
Hence, these two inputs from the actual ambipolar diffusion models guided our design of the 
experiment - we needed to sample cores and envelopes on the relevant spatial scales and 
we could (to first order) assume that the magnetic fields in the clouds were not significantly 
twisted but mainly ordered. 

We chose four clouds for our experiment - two in the Taurus molecular cloud complex 
(distance ~ 150 pc) and two in the Perseus molecular cloud complex (distance ~ 300 pc). 
A 0.1 pc core radius would be ~ l'(2') and a 2 pc cloud (or envelope) radius would be 
~ 20'(40') at the 300 (150) pc distance. At the OH line frequency, the primary-beam radius 
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of the Arecibo telescope is ~ 1.5', so this is well matched for measurement of core properties 
at the densities sampled by OH. The Green Bank telescope (GBT) beam radius is ~ 3.9'; 
by pointing the GBT at positions 6' from the Arecibo pointing position, the GBT beams 
would exclude the molecular core and sample (at half-power response) the radius range 2.1' 
to 9.9', or ~ 0.2(0.1) pc to ~ 0.9(0.45) pc at the 300 (150) pc distance. Figure 1 displays 
these beam sizes and positions. This choice of beams was set by the specifications of the 
ambipolar diffusion models cited above. The sampling of the envelope had to be sufficiently 
far from the core to obtain a significantly different result from the core result, but well within 
the outer boundary of the unevolved clouds; the GBT beamsize was ideally suited to these 
objectives. We could then "synthesize" a toroidal or ring beam to sample the envelopes by 
appropriately combining the observations from the four GBT beams in order to produce 
exactly the envelope sampling called for by the ambipolar diffusion models. 

The Arecibo OH survey of dark clo uds produced detections of Blos in 9 of 33 cores 



observed; earlier, iGoodman et al.l (119891 ) obtained a detection with Arecibo toward the Bl 
core. The Arecibo observations have therefore selected 10 cores for which there are measure- 
ments (not just limits) for Nqh/Blos (oc M/$) with a 3' beam, well matched to the core 
diameters of ~ 2'. Additional observations that are needed are measurements of Nqh/Blos 
in the envelopes of these cores to test the prediction that ambipolar diffusion increases the 
mass but not (at least very much initially) the field in cloud cores. Zeeman observations 
in molecular envelopes have generally not been attempted before, since line strengths are 
weaker and Zeeman detections are difficult to obtain. 

In order to measure M/$ in the envelopes surrounding the Arecibo cores, for this 
experiment we made a four-point Zeeman map with the GBT at positions 6' from each 
of four cores with detected B^qs (see figure 1), dividing the observing time about equally 
between the four positions. Combining the four envelope results for each cloud can then 
give results that would be obtained with the required toroidal or ring telescope beam; of 
course, sensitivity around the ring is not uniform due to four telescope beams that make up 
the ring. This sampling of the envelope regions is exactly what we require to measure M/$ 
in the envelopes, excluding the cores. Therefore, we have B^qs and Nqh separately for the 
core and envelope material. The ratios [NoH/BLos]core and [Nqh / BLos]enveiope will then 
be available for the clouds. The strong field/ambipolar diffusion theory for core formation 

requires [M / ^ core /[M/^ envelope = [NoH / BLOs]core/[NoH / BLOs]envelope > 1 (that is, that 

M/$ increases from envelope to core). So these observations could result in the ambipolar 
diffusion theory being proved wrong if the observations find this core/envelope M/$ ratio to 
be less than one. It could also result in the driven turbulence, ideal MHD simulations being 
proved wrong, if this experiment finds this core/envelope M/$ ratio to be greater than 1. It 
is of course impossible to prove that a theory is correct, but this experiment has the potential 
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to rule out one of the two competing extreme-case theories of star formation. Such a result 
would have a profound effect on further theoretical work on star formation, and would be a 
major advance in understanding the star formation process. 

We define 

^^■'' core I ^^ core /, \ 



11 = 



^''^ envelope I ^envelope 

The mass (the OH lines are optically thin) is given by 

M(xIAVA (2) 

and the magnetic flux by 

$ oc {Blos/cosO) a. (3) 

Here / is the peak intensity of the spectral line, AV is the FWHM line width, A is 
the area of sky sampled, A = Tcr"^ for the core and A = Airr"^ for the envelope, with r^ 
and Tg being the radii of the telescope beams used respectively for the core and envelope 
observations, Blos is the magnitude of the line-of-sight component of the magnetic field, 
and 6 is the angle between the line of sight and B. Since A is the same for measuring M 
and $ for the core, it divides out of the numerator of TZ; similarly in the denominator for 
the envelope. The factor cos9 also divides out. 

Theoretical calculations for ambipolar diffusion generally do not predict TZ as we have 
defined it, but rather M/ $ in the central flux tube as a function of the central density, which 



increases with time (e.g.. lCiolek fc Mouschoviaslll994l . figure 2e). For a disk morphology with 



the field along the line of sight, the central flux tube would include contributions from the 
foreground and background parts of the envelope, but these would be small and the central 
result would be dominated by the core M/$. In any case, our measurement of {M/^)core 
corresponds to what the models give. Hence, for a quantitative comparison with those 
predictions, we would need M/$ through the center of the initial cloud before ambipolar 
diffusion had created the core. Before ambipolar diffusion has acted, the initial central M/$ 
is equal to the M/$ of the initial cloud as a whole. We can infer that original central M/$ 
from our observations. The total mass and flux of the cloud do not change as the cloud 
evolves due to ambipolar diffusion, so the present {M/^)core+enveiope, obtained by adding 
together the Arecibo result for the core and the four results for the envelope, gives the initial 
central mass-to-fiux ratio (M/$)co used as the starting point of the ambipolar diffusion 
calculations. We therefore define 

^/ ^^J-corej ^core /,n 

^^■Lcore+envelope I ^^ core+envelope 
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7V p^ (M/$)c/(M/$)co, the quantity predicted by ambipolar diffusion models (e.g. JCiolek fc Mouschovias 



1994I ). The expressions in equations [2] and [3] for M and $ are used in equation |U but now 
the areas A do not divide out in the denominator and exphcit values must be used. The 
factor cos6 still divides out, however. 

Measuring relative values of M /$ {TZ or TZ') will eliminate or at least greatly reduce the 
uncertainties inherent in absolute measurements of M/$. We seek only a change in M/$ 
from envelope to core, not the absolute values themselves. This avoids all of the geometrical 
correction problems in going from B^qs to the total B and from Nobs to Nb (the column 
density along B, which is what is required to properly compute M/$). We do not need to 
know the angle cos6 between the line of sight and B, because we will be making a relative 
measurement, so cos6 in equation [3] will divide out. Similarly, no geometrical correction 
for the measured column densities will be necessary. Also, by observing B^os in the core 
and envelope using the same tracer, in this case OH, the ratio OH/H that would be needed 
to convert measured OH column density to total (H2) column density in order to find the 
absolute value of M/$ is not needed. TZ and TZ' are inferred from the directly measured 
quantities /, AV, Blos, and the radii of the telescope beams used for the core and envelope 
observations. 



4. Observations 
4.1. Arecibo Observations 

The Arecibcu observations of cores used here are not new observations, but have been 



reported al ready: Bl fiGoodman et al.lll989r) . L1544 (ICrutcher fc Trolandll2000r ). B217-2 and 
L1448CO (JTroland fc Crutcheiil2008l ). Details of these observations are described in these 
papers. These core positions were selected from catalogs of molecular peaks observed in 
CO, NH3, and other tracers. The exact positions for which OH Zeeman observations were 
carried out were refined by small Arecibo OH maps to determine the peak OH line strength 
positions. The Bl observations were carried out with the old Arecibo line feed, while the 
others used the new Gregorian feed. The spectral resolution of the Bl observations was 
lower due to the older spectrometer in use at that time, but the line is well resolved so 
this does not affect the line strength and the Zeeman measurement. The Gregorian feed 
data were originally not calibrated exactly correctly. This did not affect the determination 



^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated 
by Cornell University under a cooperative agreement with the National Science Foundation. 
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of Blos, which depends on the ratio of the Stokes I and V spectra that would have been 
affected precisely the same way. But for this study, we re-calibrated the spectra using the 
best known value of the L-band feed noise diode. This resulted in about a 7% reduction in 
the line strengths from those published previously. 



The characteristics of the Arecibo telescope beam have been carefully studied by lHeiles et al 



( l2001al ). This study was of the Gregorian feed, but the essential characteristics of the line 
feed for the purposes of this paper are not significantly different. One important parameter 
is the main beam diameter at the half-power point, 2.9'. This is the size of the core region 
sampled by the Arecibo observations. Second, there is the main beam efficiency, which is 
only 0.48. This means that slightly more than half the area of sky to which the Arecibo 
observ ations were sensitive lies outside the 2.9' diameter of the main beam, mostly in side- 



lobes. iHeiles et al.l ( l2001al ) measured the telescope response out to the first sidelobe, which 
lies ~ 5' from the beam center, close to the positions where the GBT envelope observations 
were centered. This first sidelobe has 0.14 of the spatial response of the telescope, so in addi- 
tion to the response of the Arecibo telescope to the cores, there is a response to the envelope 
emission equal to 0.14/0.48 = 29% of the core response. This first sidelobe is a "ring" beam 
similar to the pseudo-ring beam formed by the four telescope beams used for the envelope 
observations (see figure 1). Hence, for two reasons, the Arecibo results for {M/^)core are 
"contaminated" by input from the envelope region. First is the fact that the Arecibo beam 
must pass through the near and far regions of the envelope along the line of sight. This is 
okay, however, for the theoretical models predict M/$ for the central fiux tube that would 
include the near and far envelope. The second is the "ring" beam produced by the first 
sidelobe. The effect of this would be to make {M/^)core look more like {M / ^) envelope than 
its "real" value, i.e., to make TZ and TZ' closer to 1 than either would be if the core only could 
be sampled. Without much more detailed information, it is impossible to deconvolve the 
envelope contribution out of the measured {M/^)core- This means that a difference between 
a measured TZ or TZ' and 1 is more statistically significant than the significance determined 
solely from the measurement uncertainties. 



4.2. Green Bank Telescope Observations 

The observations of the envelope positions were carried out with the L-band receiver 
of the NRAol 100-m Robert C. Byrd Green Bank Telescope (GBT) between Oct 2007 and 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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Apr 2008. The main beam diameter of the GBT at 1666 MHz is 7.8' and the main beam 
efficiency is 0.81; these quantities come from our own "spider" scan polarization cahbration 
observations (see below). A negligibly small fraction of the GBT response comes from the 
core region of each cloud, due to the small sidelobe intensity, the fact that the cores lie near 
the first null of the beam response, and the fact that the core ffils only a small fraction of 
the 360° surrounding the GBT beam center. Observing time at each of the 16 positions (6' 
north, south, east, and west of each of four Arecibo cores, see figure 1) was approximately 
equal, averaging about 9.5 hours of actual on-source Zeeman integrations per position. The 
NRAO Spectral Processor was used as the backend correlation spectrometer and polarimeter. 
Simultaneous observations were made in both horizontal and vertical linear polarization of 
the 1665 and 1667 MHz OH lines. The total spectrometer bandwidth for each line in each 
polarization was 156.25 kHz or about 28 km s~^, which yielded a spectral channel width of 
about 0.055 km s~^. Bandpass calibration was carried out by frequency switching between 
-15/64 and +15/64 of the total bandwidth; after appropriately combining the two halves of 
the band, this resulted in a final velocity coverage for each spectrum of about 13 km s^^. 
The initial data processing made use of IDL software written by C. Heiles and T. Robishaw, 
which provided calibrated Stokes I, V, Q, and U spectra. Linear baselines were removed 
from the Stokes I and V spectra. The polarization calibrations used "spider" scans of the 
continuum source 3C286 over several hours surrounding its transit. iHeiles et al.l (l2001bl ) 
provide details of the polarization calibration procedure. In addition to the molecular cloud 
envelope observations, we briefly observed the polarized OH maser sources W49 and W3(0H) 
in order to verify the observing setup and data processing procedures. In addition, we 
observed the previously well-studied OH Zeeman absorption line source NGC2024 (Orion 
B, W12) in order to verify that the procedures produced the same B jos, including the 
sign or direction of the magnetic field, as had bee r i obtained previous ly (ICrutcher fc Kazes 



19831 : iHeiles fc Stevens! Il986l: ICrutcher et alJll999l: iBourke et al.ll200lr ). (Note that the field 
direction reported in the ICrutcher fc Kazed (119831 ) paper was wrong; the correct direction 
is given in the latter three papers.) The GBT Zeeman results for NGC2024 agree with the 
previous results. 

Although the Heiles-Robishaw software produces a Zeeman fit to the spectra, it does so 
over the full bandwidth. The Stokes I spectra at the molecular core positions are dominated 
by a single strong component, but at the envelope positions this component is weaker and 
additional components are present that are sometimes comparable in strength with the core 
component. We therefore used our softw are, the same code used previously to fit the Arecibo 
survey data (JTroland fc Crutcherll2008l ). to fit for B^os only over those spectral channels in 
which the strong core line cor nponent dominates; these a re the results reported here. Details 
of the procedure are given by lTroland &: Crutcherl (120081 ) and briefly described in §2. II above. 
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There are several possible methods to combine the results for the four envelope positions 
surrounding each core in order to produce a measurement of Bios for each envelope. These 
include (1) taking the mean of the four envelope results for B^oS) or (2) first averaging the 
Stokes I and V spectra for the four envelope positions before fitting for B^os- The second 
method has the disadvantage that the radial velocities and line widths are not the same at 
each envelope position, so the averaged spectra would be broader than any of the individual 
spectra, which would slightly reduce the sensitivity to B^qs- We feel that instead of either 
of these two, the best analysis technique is to fit all eight spectra (the 1667 and 1665 MHz 
spectra for the four envelope positions for each cloud) simultaneously in order to obtain the 
best-fit single value for B^os and its uncertainty from the entire envelope data set. This 
preserves the information content of each of the eight spectral lines, so no broadening takes 
place that would reduce the sensitivity. It also imposes the constraint that a single value of 
Blos describes all of the data, which is the desired result for a synthesized toroidal beam 
measurement of Blos in each envelope. 



5. Results 

Results for the core and envelope observations are shown in table 1 and in figures 2-5. 
Table 1 gives the numerical values for Iqh = ^a(1665) -|- T4(1667), A\/, and B^qs- The 
values listed for the envelopes are the results from the simultaneous fit for a single B^qs 
and its mean error over both OH lines and all four envelope positions. Plotted in each panel 
of figures 2-5 is the observed antenna temperature Ta of the 1667 MHz line at each of the 
observed positions, in order to show the relative line strengths and widths at the core and 
the envelope positions. The 1665 MHz line strengths are typically about 60% of the 1667 
MHz line strengths and show the same relative strengths from position to position. In the 
upper left of each panel is Blos and its la uncertainty for that position. Although the line 
strength does not decrease uniformly in all directions from the cores, and in most cases there 
is OH emission at velocities slightly away from the velocity of the main core component, 
the line strength in the envelope is typically about 50% that of the core. Moreover, in all 
cases, the line width at the core position is only about 75% that at the envelope positions. A 
decrease in line width from envelope to core is a standard feature of molecular clouds, seen 
even more strongly in molecular tracers that sample higher densities than does OH. So both 
the increase in line strength and the decrease in line width from the GET to the Arecibo 
spectra argue that the Arecibo OH observations do sample molecular cores, and that the 
GET observations are dominated by non-core emission. 

Blos at L1544west is twice Blos at L1544core in spite of the OH column density ratio 
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being 0.6; this is the only case of a stronger field at an envelope position than in the core. 
Toward Bl there appears to be a north-south ridge of B^os with much weaker Blqs to the 
east and west of the core. At all of the remaining 13 envelope positions the signal/noise ratio 
is low, although at all 16 positions the sensitivity to B^qs is sufficient to provide significant 
results for TZ and TV . A hypothetical measured envelope Blqs = is a perfectly acceptable 
experimental result for computing TZ and TZ'\ the only requirement is that the sensitivity to 
Blos be sufficient to yield statistically meaningful results for TZ and TZ' . 

Table 2 shows TZ and TZ! and the la uncertainties for each cloud. The uncertainties in 
IZ and IZ' depend on the uncertainties for the B^os given in table 1 and the uncertainties in 
Iqh and ^V . The uncertainty in the A^ is about 0.02 km s~^ in each case. The nominal 
uncertainty in Iqh from the channel-to-channel noise in the spectra is only about 0.01 K. 
The uncertainties in both Iqh and AV are too small to contribute to the uncertainties in TZ 
and TZ' . However, the absolute calibration of the line strengths is uncertain by about 10%. 
Although this systematic uncertainty would divide out in TZ and TZ' if the same telescope 
had been used for all observations, it may be that there is a calibration difference between 
the Arecibo telescope and the GBT. We have therefore used an uncertainty of 10% in each 
of the core and envelope values for Iqh- Even so, this uncertainty is insignificant in the error 
budgets for TZ and TZ' . Those error budgets are dominated by the uncertainty in the B^os 
of each envelope, since the signal/noise ratio for those quantities is by far the lowest of any 
of the four measured quantities that go into TZ and TZ' . We assume that the PDFs of these 
four measured quantities are Gaussian normal distributed (see §2.11) . 

We obtained the results for TZ and TZ' and their la uncertainties with two different 
methods. First, we compute the uncertainties in TZ and TZ' by normal error propagation. 
This procedure is justified by the fact that the uncertainties in Blqs are Gaussian normal (or 
very nearly so), and the possible systematic uncertainty in the N contribute insignificantly 
in comparison with the uncertainties in the Blqs- However, due to the low signal-to-noise 
ratio in the B^os (especially in the envelope field strengths), the assumption of Gaussian 
normal errors in TZ and TZ' may not be strictly correct. We have therefore computed the 
PDF for TZ and TZ' for each of the four clouds by Monte Carlo simulations; for each Monte 
Carlo simulation, 10^ trials were used. The results for TZ are shown in figure 6; results for 
TZ' are extremely similar. Note that TZ and TZ' may be positive or negative. Over plotted on 
the Monte Carlo PDFs are the Gaussian normal error curves with the parameters given by 
the error propagation method. In all cases, the means and standard deviations from the two 
methods are essentially identical. However, the Monte Carlo PDF generally has a stronger 
tail at the high values of TZ and TZ' ^ which slightly increases the probability that TZ (and TZ') 
are greater than one. These probabilities are essentially identical for TZ and TZ' , so we list 
only one probability for each cloud in table 2. Each of the four clouds is unlikely to have TZ 
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and 7?.' > 1, with probabilities of being greater than one varying from less than 1% to about 
10% (table 2). The probability that all four of our clouds have 7^' > 1 is just the product 
of the probabilities in table 2, or 3 x 10~^, a highly significant result. Our experiment is 
therefore in contradiction with the hypothesis that these four cores were formed by ambipolar 
diffuson. 

There are several possible biases that could affect 71 and 71'. First, as discussed above, 
the Arecibo telescope beam pattern biases the observed 7Z and 7?.' to be closer to 1 than 
the "real" values. Since the observed 7^ and 7^' < 1, these observed values are biased to be 
higher than the actual values. In addition, there is a factor that may systematically raise 
or lower 7^ and 7^' - the possible curvature of the magnetic field lines. Field lines will be 
drawn into an "hourglass" morphology as a core forms. Because the four cores out of 34 that 
we have observed in this experiment are among the few with strong Blqsi it is likely that 
B points approximately along the line of sight for these four cores. The field line through 
the cloud center would then point along the line of sight, so Blqs would equal the total B. 
However, other field lines passing through a core would be curved away from the line of sight 
in the near and far side of the core, so Blos would be less than the total B. In the envelope 
a similar effect would hold, but the curving of magnetic fi eld lines in the "ho urglass" would 



be less extreme in this lower density outer region. See I G alii &: Shd (Il993l ) figure 5b for 
an example of the morphology of magnetic field lines; for a typical distance of 150 pc, the 
radius of the Arecibo beam would be 2 x 10^^ cm. Hence, one might generally expect Blos 
to underestimate -B by a smaller factor in the envelope than in the core of a cloud, which 
would bias 7^ and 7?.' to be too high - the same sense as the Arecibo sidelobe bias. However, 
exactly which direction the field curvature bias would go and by how much would depend 
on the detailed structure of the magnetic field morphology. But in any case, the bias must 
be insignificant. Even if the mean field directions over the GBT and Arecibo beams differed 
by 30° due to an hourglass morphology, a very unrealistically large value, the change in the 
measured Blqs would be by the factor cos30° = 0.87, and 7^ and 7^' would be biased high 
by the factor 1/0.87. The resulting 7^ and 7^' would be changed by only about la by such an 
unrealistically large net angle difference. (Of course, if the magnetic fields are not regular, 
this angle could be much larger, but for our comparison with the "idealized" ambipolar 
diffusion model this is not relevant.) Finally, since we assume that N{OH) oc the mass M, if 
[OH/H] varies between the envelope and core, 7?. and 7Z' would be affected. It seems unhkely 
that there is a significant systematic variation in [OH/H] between the envelopes and cores 
that we have observed. There is only a 50% difference in the column density of OH between 
core and envelo pe, so the physical regions being sampled are not that different. Moreover, 
Crutcherl (119791 ) found no evidence for a variation in [OH/H] up to n{H2) ~ 2 x 10^ cm~'^. 



which includes the range our observations sample. So our measured 7Z and 7Z' can be slightly 
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biased, probably overall in the direction of being closer to 1 than the "real" values, but not 
by sufficiently large amounts to change the significant result of this experiment. 



6. Discussion 

We have found that M/$ decreases significantly from envelope to core, or from the initial 
central value to the present evolved central value, for the clouds we have studied. How does 
this observational result compare with the prediction of the strong field, ambipolar diffusion 
driven theory and the weak field, turbulence driven theory? 

We first consider the ambipolar diffusion theory. We compare our results with what 
we called above idealized models, which include only gravity, regular magnetic fields, and 
thermal pressure. The non-thermal motions and irregular structures of observed clouds, 
which themselves are embedded in a complex interstellar medium, are not directly considered 
in these idealized models. Perhaps for this reason, the ambipolar diffusion models are not 
compatible with the structure we observe in both column density and magnetic field strength. 
Figures 2-5 show clearly that the column densities vary considerably around the core positions 
of our clouds, rather than the uniform result that would be expected from a idealized model 
with B oriented closely along the line of sight. Also, the observed Blqs in two envelopes 
clearly show structure, with B^os at the west envelope position of L1544 being about twice 
as strong as toward the core, and much stronger than toward the other three envelope 
positions. Bl shows a north-south "ridge" of Blqs, with weaker Blqs toward the east and 
west envelope positions. These observations cannot be reproduced by an ambipolar diffusion 
idealized model. 

The strong magnetic field/ambipolar diffusion theory requires that M/$ increase in the 
core as evolution proceeds; after all, this increase is ambipolar diffusion, the heart of this 
theory. Hence, this theory predicts 7?.' > 1. The amount by which 7?.' exceeds 1 would be 
dependent o n the specific parameters for a model cloud - mainly the assumed initial M/$. 



For example, ICiolek &: MouschoviasI (119941 ) discussed a model cloud broadly consistent with 
the parameters of the clouds we observe. Although their model is in terms of dimensionless 
parameters, they state that the original unevolved cloud is consistent with a temperature of 
10 K, a central density of 2.6 x 10^ cm~^, a central magnetic field of 35 fxG, a radius of 4.3 
pc, and a total mass of 98 Mq. At the time when the core becomes magnetically critical, 
the central density has increased by a factor ~ 37 while the central B has increased by less 
than 1.7. The initial M/$ is subcritical with the central M/$ = 0.256 of critical; for this 
model, the predicted TZ' ~ 4. However, our results for all four observed clouds is 7^' << 4 
with high degrees of significance. 
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There have been ambipolar diffusion models calculated specifically for t wo of the clouds 



i n our sample, for comparison with observational data available at that time. ICrutcher et al. 



(119941 ) discussed a model specifically for Bl; it had a core mass of 13 Mq, an envelope 
mass of 600 Mq, an envelope radius of 2.9 pc, an initial central i? = 43 iiG, and an initial 
central M/$ = 0.42 of critical. The model assumed that the cloud was a disk whose minor 
axis was at an angle 9 = 70° to the line of sight; all observed properties of Bl available 
at that time were given accurately by the model. The prediction of this model would be 
TZ' = 1/ 0.42 = 2.4, a factor of 5 l arger than our result (table 2) for Bl. Moreover, now 



that the iTroland fc Crutcherl (120081 ) survey of dark cloud cores has shown that Bl has the 
greatest B^os of any core with a detected B^osi it seems more likely that B is nearly along 
the line of sight, and that the true central total B is close to the observed B^os = 27 /xG and 
not the model result i? = 85 /iG, which implied B^os = 85cos70° = 29 fiG. When Bl was 
among a very small number of dark clouds with sensitive OH Zeeman observations, it was 
not unreasonable to hypothesize that its field lay nearly in the plane of the sky. However, 
other clouds similar to Bl with similar total field strengths sho uld have B nearly along the 



line of sight, yielding B^qs ^ 85 /xG; these are not found in the lTroland fc Crutcherl (120081 ) 
survey results. 



Ciolek &: Basul (120001 ) computed a model for L1544. The model had 30 Mq within a 
radius of 0.45 pc, with additional mass that does not participate in the evolution in the 
envelope at larger radius (not clearly specified, but apparently about 2.5 pc based on their 
figure Ic). They assumed 6 ~ 74°, again a very large angle between B and the line of 
sight, which was necessary in order to have the required large central B agree with the small 
observed value of Blos- The initial M/$ was 0.8 of critical, or closer to critical than other 
ambipolar diffusion models discussed above. This would imply TZ' ^ 1.25, which differs from 
our measurement 7?.' = 0.46 ± 0.43, although not by a highly significant amount. However, 
as for Bl, the large value for 6 they had to assume in order to make the field strength of the 
model agree with the observation of B^qs seems unreasonably large. Not all B can lie near 
the plane of the sky! 

The ambipolar diffusion model results are all 7^' > 1, with the actual value depending 
on the initial assumed M/<l>. But even if all clouds start only very slightly subcritical, which 
would in itself minimize the importance of ambipolar diffusion in cloud evolution, our results 
are not consistent with the ambipolar diffusion requirement. 

A possible way out of this conclusion might be to hypothesize that the magnetic fields 
in our envelope regions are not the regular ones of the idealized models, but rather that 
these magnetic fields twist and indeed reverse direction. By fitting for a single value of 
Blos for each envelope, we obtain a smaller value for Blos (and hence a smaller TV) in 



-19- 



the entire envelope surrounding the core than we would have gotten if we had (for example) 
just taken the largest value of B^os at any of the four positions. If the field at the GBT 
positions twisted significantly, then perhaps the largest absolute value of B^os would be the 
appropriate one to use in calculating the magnetic flux in the envelope region. At the other 
three positions surrounding a core, Blqs could be smaller due to B being twisted into the 
plane of the sky or even to pointing in the opposite direction from the field in the core. 
However, our fit for a single envelope value of B^qs is the proper one for testing at least 
the idealized models of this theory. This single value gives (approximately) the result that 
would be obtained with a ring telescope beam as it would sample the envelope regions. But 
more importantly, the GBT beams are centered only 6' from the center of each core, or at 
a radius of 0.26 pc for a typical cloud distance of 150 pc. This radius is about an order 
of magnitude smaller than the cloud radii in the three idealized ambipolar diffusion models 
discussed above, where the models require a smooth and not tangled magnetic field. The 
GBT positions are not far enough away from the cores that they could be sampling unrelated 
flux tubes that may have nothing to do with the envelopes of the cores. Moreover, such a 
twisted (indeed reversed) field explanation would have to hold for all four of our cores, which 
seems statistically unlikely. Finally, the good agreement in radial velocities of the core and 
envelope OH spectral lines (figures 2-5) suggests that the GBT observations do sample the 
envelope regions of the cores. 

What about predictions of initially weak magnetic field, turbulent theories? Such simu- 
lations are of course also idealized, with artificial boundary conditions, simple turbulence 
algorithms without a physical mechanism for generating the turbulence, and neglect of 
some physics. None theless, our comparison is with state-of-the-art simulations. Recently, 



Lunttila et al.l (120081 ) have computed synthetic Zeeman profiles for their super- Alfvenic (weak 
magnetic field) simulations of molecular clouds and cores formed by turbulence in an ini- 
tially uniform interstellar medium. The mean field strength in their simulations is quite low; 
Blos ~ 2.1 /iG along the direction of the original uniform B. They performed numerically 
the same experiment for which we report here the observational results. That is, they "ob- 
served" the synthetic Stokes I and V spectra at each of their core positions with a 3' beam 
and at four envelope positions around each core with 8' beams pointed 6' north, south, east, 
and west of the core positions. They computed N{OH) from their Stokes I profiles and Blqs 
by fitting their synthetic Stokes dl /dp to their V profiles (the same technique used obser- 
vationally). They analyzed 36 (out of 139 total) cores with B^qs > 10 /iG, corresponding 
roughly to field strengths in the four cores we have observed. For these 36 cores, the mean 
N{H2) = 6.6 X 10^^ cm~^, which agrees reasonably well with the mean N{H2) = 4.3 x 10^^ 
cm~^ for our four cores. The mean line width of the 36 cores was AV ~ 1 km s^^, which also 
agrees with our four observed cores. For these cores they found that TZ had a large scatter. 
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0.08 < 7^ < 1.6, with the PDF favoring values less than 1. Hence, the weak field, turbulent 
calculation of the formation of molecular clouds and cores seems to agree in physical prop- 
erties with the four cores we have observed. Given that our experiment included only four 
clouds, the range in TZ corresponds well with the observed range. 

We suggest that other simulations, such as those with stronger magnetic fields that 
include both ambipolar diffusion a nd turbulence, s hould be compared with our experimental 



results in the same manner as did iLunttila et al.l (120081 ). 



7. Conclusion 

Previous Zeeman studies of magnetic fields in molecular clouds have not been definitive 
in testing the two extreme-case models of star formation. The mean mass-to-fiux ratios 
M/$ found from these statistical studies were slightly supercritical - consistent with either 
theory. Detailed ambipolar diffusion models for two clouds found excellent agreement with 
the observations, although both required the field to be nearly in the plane of the sky in 
order not to produce line-of-sight fields much stronger than observed. Uncertainties in the 
angle between B and the line of sight and in the total hydrogen column density are inherent 
in measuring with Zeeman observations. In order to mitigate these uncertainties, we have 
measured the ratio of M/$ between the envelopes and cores of four molecular clouds in order 
to test ambipolar diffusion (strong magnetic fields) versus turbulence (weak magnetic fields) 
driven star formation theory. The theory of star formation that hypothesizes clouds initially 
supported by strong magnetic fields, with evolution and core formation being driven by 
ambipolar diffusion, predicts that the central M/<l> must increase as ambipolar diffusion acts. 
Idealized models predict that the increase in M/<l> up to the point when the core becomes 
supercritical and gravitational collapse proceeds is approximately equal to the inverse of the 
amount by which the original cloud was subcritical; that is, TV > 1. The probability that 
all four of our clouds have 7^' > 1 is 3 x 10^^, a highly significant result. On the other 
hand, simulations which form clouds and cores by turbulence acting in a weak magnetic 
field environment preferentially yield a M/$ ratio between core and envelope TZ < 1, in 
agreement with our results. 

Telescope availability limitations allowed only four clouds to be observed; unfortunately 
the extremely large amount of telescope time required precludes expanding this experi- 
ment beyond four clouds for the foreseeable future. The theoretical predictions of TZ and 
TZ' are based on idealized ambipolar diffusion models and idealized turbulence simulations. 
Nonetheless, the clear conclusion from our experiment is that at least for these four clouds, 
the prediction of the idealized ambipolar diffusion models does not agree with our obser- 
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vational results, while the prediction of initially supercritical turbulence-driven simulations 
does. Still untested is whether simulations that include both significant magnetic fields and 
turbulence better match the data than either of the extreme cases. We suggest that all 
theorists who simulate the formation and evolution of molecular clouds and cores test their 



simulations against the results of this experiment in the manner of iLunttila et al.l (120081 ): 
that is, by calculating Stokes I and V spectra of OH from the simulations and "observing" 
Blos with our beam patterns (figure 1). 
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hours of observing time on the GET in response to proposal GBT07A-029, and the obser- 
vatory staff for help in making the observations successful. We thank Phil Perillat of the 
Arecibo Observatory for help in re-calibrating the Arecibo spectra, and we especially thank 
Carl Heiles and Tim Robishaw for allowing us to use their GET software and for helping us 
understand how to make it work. Finally, we thank two conscientious referees who provided 
valuable suggestions to refine the scientific and statistical arguments in this paper. This 
research was partially supported by NSF grants AST 0307642 and 0606822, and by NRAO 
grant GSSP07-0007. 

Facilities Arecibo, GET 



REFERENCES 

Eourke, T. L., Myers, P. C., Robinson, G., & Hyland, A. R. 2001, ApJ, 554, 916 

Ciolek, G. E. & Mouschovias, T. Ch. 1994, ApJ, 425, 142 

Ciolek, G. E. & Easu, S. 2000, ApJ, 529, 925 

Crutcher, R. M. 1979, ApJ, 234, 881 

Crutcher, R. M. & Kazes 1983, A&A, 125, 23 

Crutcher, R. M., Troland, T. H., Goodman, A. A., Heiles, C, Kazes, I., & Myers, P. C. 
1993, ApJ, 407, 175 

Crutcher, R. M., Mouschovias, T. Ch., Troland, T. H., & Ciolek, G. E. 1994, ApJ, 427, 839 

Crutcher, R. M., Roberts, D. A., Troland, T. H., & Goss, W. M. 1999, ApJ, 515, 275 

Crutcher, R. M. 1999, ApJ, 520, 706 



-22- 

Crutcher, R. M. & Troland, T. H. 2000, ApJ, 537, 139 

Elmegreen, B. G. 2000, ApJ, 530, 277 

Elmegreen, B.G. & Scalo, J. 2004, ARAAp, 42, 211 

Gain, D. & Shu, F. H. 1993, ApJ, 417, 220 

Goodman, A. A., Crutcher, R. M., Heiles, C, Myers, R G., & Troland, T. H. 1989, ApJ, 
338, L61 

Heiles, G. & Stevens, M. 1986, ApJ, 301, 331 

Heiles, G. & 10 coauthors 2001a, PASP, 113, 1247 

Heiles, G. & 9 coauthors 2001b, PASP, 113, 1274 

Heiles, G. & Grutcher, R. M, in Cosmic Magnetic Fields, ed. R. Wielebinski & R. Beck, R. 
(Springer: Berlin, Heidelberg), 2005 (Springer: Berlin, Heidelberg), 2005 

Kudoh, T. & Basu, S. 2008, ApJ, 679, L97 

Lunttila, T., Padoan, P., Juivela, M., & Nordlund, A. 2008, ApJL in press (astro-ph 
0806.3854) 

MacLow M.-M. & Klessen R. S. 2004, Rev. Mod. Phys., 76, 125 

McKee, G. F. & Ostriker, E. G. 2007, ARAAp, 45, 565 

Mouschovias, T. Gh. & Spitzer, L., Jr. 1976, ApJ, 210, 326 

Mouschovias, T. Gh. 1991, in The Physics of Star Formation and Early Stellar Evolution, 
NATO ASIG Proc. 342, ed. G. J. Lada & N. D. Kylafis., p. 61 

Mouschovias, T. Gh. & Giolek, G. E. 1999, in The Origin of Stars and Planetary Systems, 
ed. G. J. Lada & N. D. Kylafis (Kluwer: Dordrecht), p. 305 

Nakano, T. & Nakamura, T. 1978, PASJ, 30, 681 

Nakamura, F. & Li, Z.-Y. 2005, ApJ, 631, 411 

Nakamura, F. & Li, Z.-Y. 2008, ApJ, accepted 

Padoan, P. & Nordlund, A 1999, ApJ, 526, 279 

Padoan, P., Jimenez, R., Juvela, M., & Nordlund, A 2004, ApJ, 604, L49 



-23- 



Shu, F., Allen, A., Sliang, H., Ostriker, E. C, & Li, Z.-Y. 1999, in The Origin of Stars and 
Planetary Systems, ed. Charles J. Lada & Nikolaos D. Kylafis, (Kluwer: Dordrecht), 
p. 193 

Troland, T.H. & Crutcher, R.M. 2008, ApJ, 680, 457 

Vazquez-Semadeni, E., Kim, J., Shadmehri, M., & Ballesteros-Paredes, J. 2005, ApJ, 618, 
334 



This preprint was prepared with the AAS I^TJtjX macros v5.2. 



-24- 



Table 1. Observational Results 



Cloud 


Ioh{K) 


^V{km/s) 


Blos{i^G) 


L1448CO(env) 


0.63 


1.16 


-0±5 


L1448CO(core) 


1.30 


0.93 


-26 ±4 


B217-2(env) 


0.59 


0.79 


+2 ±4 


B217-2(core) 


1.28 


0.47 


+14 ±4 


L1544(env) 


0.96 


0.67 


+2 ±3 


L1544(core) 


2.43 


0.48 


+11 ±2 


Bl(env) 


1.21 


1.32 


-8 + 3 


Bl(core) 


1.93 


1.14 


-27 + 4 



Table 2. Relative Mass/Flux 



Cloud 


n 


W 


Probability 7^ or 7^' > 1 


L1448CO 


0.02 + 0.36 


0.07 + 0.34 


0.005 


B217-2 


0.15 + 0.43 


0.19 + 0.41 


0.05 


L1544 


0.42 + 0.46 


0.46 + 0.43 


0.11 


Bl 


0.41 + 0.20 


0.44 + 0.19 


0.010 
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Fig. 1. — The Arecibo telescope primary beam (small circle centered at 0,0) and the four 
GBT telescope primary beams (large circles centered 6' north, south, east, and west of 0,0. 
The dotted circles show the first sidelobe of the Arecibo telescope beam. All circles are at 
the half-power points. 
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Fig. 2. — OH 1667 MHz spectra toward the core of L1448CO obtained with the Arecibo 
telescope (center panel) and toward each of the envelope positions 6' north, south, east, and 
west of the core, obtained with the GET. In the upper left of each panel is the inferred B^os 
and its la uncertainty at that position. A negative Blqs nieans the magnetic field points 
toward the observer, and vice versa for a positive B^os- 
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Fig. 3.— As in figure 2, but for B217-2. 
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Fig. 4. — As in figure 2, but for L1544. 
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Fig. 5. — As in figure 2, but for Bl. 
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Fig. 6. — Monte Carlo estimate of TZ for each cloud. The Gaussian normal error curve for 
the error propagation results given in table 2 are also shown. 



